Center for Turbulence Research 
Annual Research Briefs 2003 


251 


The effect of the density ratio on the nonlinear 
dynamics of the unstable fluid interface 

By S.I. Abarzhi 


1. Motivation and objectives 

When a light fluid accelerates a heavy fluid, the misalignment of the pressure and 
density gradients gives rise to the instability of the interface, and produces eventually 
the turbulent mixing of the fluids (Rayleigh 1892; Davies & Taylor 1950; Richtmyer 1960; 
Meshkov 1969). This phenomenon is called the Rayleigh-Taylor instability (RTI) if the 
acceleration is sustained, and the Richtmyer-Meshkov instability (RMI) if the acceleration 
is driven by a shock or if it is impulsive. The RT/RM turbulent mixing is of extreme 
importance in astrophysics, inertial confinement fusion, and many other applications 
(Sharp 1984). To obtain a reliable description of the mixing process, the evolution of a 
large-scale coherent structure, the dynamics of small-scale structures, and the cascades 
of energy should be understood. 

The large-scale coherent structure is an array of bubbles and spikes periodic in the plane 
normal to the direction of gravity or the initial shock (Rayleigh 1892; Davies & Taylor 
1950; Meshkov 1969; Schneider et al. 1998). It appears in the nonlinear regime of RTI and 
RMI and has a spatial period determined by the mode of fastest growth (Chandrasekhar 
1961). The light (heavy) fluid penetrates the heavy (light) fluid in bubbles (spikes). The 
density ratio is a determining factor of the instability dynamics (Sharp 1984; Schneider 
et al. 1998). Singular aspects of the interface evolution (such as the generation of vorticity 
and secondary instabilities, resulting in the direct and inverse cascades of energy, a finite 
contrast of the fluid densities and the non-linearity of the dynamics) cause theoretical 
difficulties and preclude elementary methods of solution (Dalziel et al. 1999; He et al. 
1999; Gardner et al. 1988; Jacobs & Sheeley 1996; Kucherenko et al. 2000; Holmes et al. 
1999; Volkov et al. 2001). 

For fluids with highly contrasting densities (fluid- vacuum) , the effect of singularities 
on the interplay of harmonics and on the nonlinear motion in RTI/RMI has been studied 
intensively over the decades (Layzer 1955; Garabedian 1957; Tanveer 1993). A new ap- 
proach based on group theory has been developed recently by Abarzhi (1998, 2002). The 
asymptotic theories of Layzer (1955); Garabedian (1957); Abarzhi (1998, 2002) agreed 
with experiments and simulations. For fluids with a finite density contrast, the influ- 
ence of singularities on the cascades of energy and the large-scale coherent dynamics 
in RTI/RMI has yet to be elucidated (Baker et al. 1982; Moore 1979; Hou et al. 1997; 
Cowley et al. 1999; Matsuoka et al. 2003). The empiric models proposed by Sharp (1984) 
and Oron et al. (2001) could not explain observations completely and were a subject 
for controversy. These models disregarded the conservation of mass and introduced ad- 
justable parameters to balance drag, buoyancy, and inertia in the flow, (Dimonte 2000). 
In a recent attempt of Goncharov (2002) to reproduce the results of the drag model 
of Oron et al. (2001) in a single-mode approximation, a complete set of the boundary 
conditions were not satisfied (see below), and the conservation laws were thus violated. 

Here we report multiple harmonic theoretical solutions for a complete system of conser- 
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vation laws, which describe the large-scale coherent dynamics in RTI and RMI for fluids 
with a finite density ratio in the general three-dimensional case. The analysis yields new 
properties of the bubble front dynamics. In either RTI or RMI, the obtained dependen- 
cies of the bubble velocity and curvature on the density ratio differ qualitatively and 
quantitatively from those suggested by the models of Sharp (1984), Oron et al. (2001), 
and Goncharov (2002). We show explicitly that these models violate the conservation 
laws. For the first time, our theory reveals an important qualitative distinction between 
the dynamics of the RT and RM bubbles. Asymptotically, the RT bubble is curved, and 
its curvature has a strong dependence on the density ratio, while the RM bubble flattens 
independently of the density ratio. The velocity of the RT bubble depends has a power- 
law dependence on the bubble curvature with exponent 3/2 and a universal coefficient 
independent of the density ratio, while the RM bubble decelerates. The bubble curvature 
and velocity depend mutually on one another, as do the differences between the RM and 
RT cases for these quantities. Our theory explains existing data, formulates the univer- 
sal properties of the RT and RM nonlinear dynamics, and identifies sensitive diagnostic 
parameter for future observations. 


2. Governing equations 

Let t be time, ( x , y, z ) be the Cartesian coordinates, and 9(x, y , z, t) be the scalar 
function with 9 = 0 at the fluid interface. Locally, 6 = z*(x,y,t) — 2 where z*(x,y,t) 
is the position of the fluid interface. The fluid density and velocity have the form p = 
phH{—6) + piH{9), and v = VhM(—9) + viH(9 ), where H is the Heaviside step-function, 
and ph.(i) and Vhg) are the density and velocity of the heavy (light) fluid located in the 
region 9 < 0 (9 > 0). For incompressible fluids, V • v = 0, and the values of p /,p) are 
independent of the coordinates and time. The equation of continuity is reduced then to 

(0 + v h -V9)p h \ e=o = (0 + « r V0)«| e=o (2.1) 

where the dot indicates a partial time-derivative. If there is no mass flux across the 
moving interface, the normal component of velocity is continuous at the interface and 

v h -V9\ 9=0 = Vl -V9\ e=0 = -e\ e=0 . (2.2) 

With neglected terms for viscous stress and surface tension, the momentum equation 
is transformed into the conditions 


Ph ( Vh + (v h ■ V)v/J + g) + S7p h e<0 = 0 , 

(2.3) 

pi (in + (vi • V)n ; +g) + X7 Pl 0>Q = 0, 

(2.4) 

Ph ~Pl \ g= o = 

(2.5) 


where Ph(i) is the pressure of the heavy (light) fluid, and g is the gravity directed from 
the heavy fluid to the light fluid with \g\ = g. There are mass sources in the flow, and 
the boundary conditions at the infinity close the set of the governing equations 

v h\e=-oo = ^le= + oo = °- ( 2 - 6 ) 

The spatial period, the time-scale, and the symmetry of the motion in (2. 1-2.6) are 
determined from the initial conditions. We choose the spatial period A in a vicinity of 
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the wavelength of the mode of fastest growth, A ~ A max , where A max is set by surface 
tension and viscosity (Chandrasekhar 1961). For RTI, the time-scale is trt ~ \J A / Ag , 
where A = (ph — pi)/(ph + Pi) is the Atwood number with 0 < A < 1. For RMI, 
g = 0 and trm ~ A/fo, where vo is the initial velocity value. Based on the experimental 
observations, we separate scales and divide the fluid interface into active and passive 
regions, similarly to Aref & Tryggvason (1989). In the active regions (scales <C A) the 
vorticity is intensive, while the passive regions (scales ~ A) are simply advected. In order 
to be stable under modulations with scales A, the large-scale coherent motion must 
be invariant under one of symmorphic groups with inversion in the plane (Abarzhi 1998, 
2002). 

To describe the dynamics of the nonlinear bubble, we reduce (2. 1-2.6) to a local dy- 
namical system. All calculations are performed in the frame of reference moving with 
velocity in the ^-direction, where v(t) is the velocity at the bubble tip in the laboratory 
frame of references. For the large-scale coherent motion, Vh(i) = and in the case 

of a 3D flow with hexagonal symmetry 


oo 3 

&h = d) m (t) (z + (exp(— mkz)/3mk) ^ cos (mfejr)) + cross terms + fh(t), (2.7) 

m= 1 i= 1 

oo 3 

^ l» TO (t) (— z(exp(— mkz) / 3mk ) cos (mkjr)) + cross terms + fi(t), (2.8) 

m— 1 i=l 

where fc; are the vectors of the reciprocal lattice, r = ( x,y ), and are time- 
dependent functions. For x ss 0 and y ss 0 the interface can be expanded as a power 
series, z* = Y1n= i C N{t){x 2N + y 2N ) + cross terms, where £i (t) is the principal curvature 
at the bubble tip, and £i(t) < 0. 

Substituting these expressions in (2. 1-2. 6), taking the first integral of (2.5), and re- 
expanding (2. 1-2.6) x, y « 0, we derive a dynamical system of ordinary differential equa- 
tions for the variables Cat and the moments M n (t) = '£2™ =1 $ m {t){km) n + cross terms 
and M n (t) = ^m{t){km ) n Across terms, where n is an integer. The moments are 

correlations functions by their physical meaning. For N = 1, the conditions in (2. 2-2. 6) 
take respectively the form 


Cl = 2Ci Mr + M 2 /4 = 2CiMr - M 2 /4, (2.9) 

(Mi / 4 + Ci M 0 - M\l 8 - Ci g)p h = (A/4 - Ci^o - Mr/8 - Ci g)pi, (2-10) 

M 0 (t) = -M 0 (t) = -v(t) . (2.11) 

The local dynamical system (2.11) describes the dynamics of the bubble in a vicinity 
of its tip as long as the spatial period A of the coherent structure is invariable. The 
presentation in terms of moments M n and M n allows one to perform a multiple harmonic 
analysis and find the regular asymptotic solutions with a desired accuracy. 


3. Regular asymptotic solutions 

Retaining only the first order amplitudes in the expressions for the moments, we derive 
a nonlinear solution of the Layzer-type, which conserves mass, momentum, and has no 
mass sources. In RTI, for t/rnT ^ 1 the curvature and velocity of the Layzer-type 
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bubble are Ci = Cl — —Ak/8, and v = Vl.rt = \J Ag/k, in agreement with the empiric 
approach of Sharp (1984). In RMI, for t/r 1, the curvature and velocity are Ci = 
Cl = —Ak/8, and v = Vl.rm = (2 — A 2 )/Akt. These single-mode solutions however 
do not satisfy (2.2) and (2.9) and permits mass flux across the interface. To avoid this 
difficulty Goncharov (2002) has violated the boundary conditions in (2.6) and (2.11), and 
introduced an artificial time-dependent mass flux of the light fluid in the flow. For the 
solutions of Goncharov (2002), Ci = Cn = —k/8 , and v = vd.rt = \Z^Ag/(l + A)k in 
RTI, and Ci = Cd = —k/8, and v = Vd.rm = 2/(1 + a)kt in RMI, in agreement with 
drag model of Oron et al. (2001). We conclude that Layzer-type approach (either in our 
version or in the models Oron et al. (2001) and Goncharov (2002)) does not satisfy the 
complete set of the conservation laws, and the single-mode solutions are not therefore 
physical. The reason for this difficulty lies in a non-local character of the non-linearity 
in (2. 1-2.6) and (2.9-2.11). 

To find regular asymptotic solutions, describing the nonlinear evolution of the bubble 
front in RTI or RMI, one should account for non-local properties of the flow that has 
singularities (Abarzhi et al. 2003). The singularities determine the interplay of harmonics 
in the global flow as well as in the local dynamics system. They transfer the fluid energy to 
smaller and larger scales and generate higher order harmonics. If the energy transports 
are not extensive, so the symmetry and the spatial period of the coherent structure 
do not change, the singularities affect the shape and velocity of the regular bubble. 
Assuming the bubble shape, parameterized by the principal curvature(s) at its tip, is 
free, we find a continuous family of regular asymptotic solutions for the local system. 
The family involves all solutions allows by the symmetry of the global flow. For the 
regular asymptotic solutions the interplay of harmonics is well captured. We perform a 
stability analysis and choose the fastest stable solution in the family as being physically 
significant. The reader is referred to Abarzhi et al. (2003) for more details. 

For the Rayleigh-Taylor instability, the bubble velocity is the function on the bubble 
curvature and density ratio, v = v(\J g /k. A, (Ci /&))• In the interval Ccr < C < 0, the 
Fourier amplitudes <I> m and <f> m decay exponentially with increase in their to, the lowest- 
order amplitudes are dominant, and higher order corrections for family solutions are 
small. For / ~ Ccr the convergence is broken, where Ccr ~ — fc/6 for A ss 1 and Ccr ~ 0 
for A ss 0. For fluids with highly contrasting densities, A ~ 1, the magnitudes of the 
Fourier amplitudes of the light fluids are much larger that those of the heavy fluid, so 
<l > rn | ~ |$ m+1 |. For fluids with similar densities, A ~ 0, this difference is insignificant, 
and |<F m | ~ |$ m |. 

The fastest solution in the family has the curvature and velocity 

Ci = C a,rt, v = va.rt (3.1) 

Explicit analytical expressions for Ca.rt and va.rt are cumbersome and not presented 
here. Remarkably, the bubble velocity and the curvature obey a universal dependence 

va.rt = y/ g/k{8\CA,RT I /k) 3/2 - (3.2) 

In the limiting case of fluids with highly contrasting densities, A ss 1, the solution (3.1) 
takes the form 

Ca.rt ~ — (&/8)(l — (1 — A)/8), va.rt- v^A(1-3(1-A)/16). (3.3) 

In the other limiting case, A « 0, it can be expanded as 
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Ca.rt ~ -~(k/2)A 1 / 3 , v AiB T » (3/2) 3/: WAg/3k, (3.4) 

Stability analysis shows that solutions with ( ~ ( cr and ( ~ 0 are unstable, while 
solutions with ( ~ Ca,rt are stable. Therefore, the physically significant solution in the 
RT family is the solution (3.1). 

For the Richtmyer-Meshkov instability, the multiple harmonic regular asymptotic so- 
lutions in (2.9-2.11) with M n ~ 1/t, M n ~ 1/t and v ~ 1/f, and time-independent C/v 
are found in a similar way. The RM and RT families have a number of common prop- 
erties such as the dominance of the lowest order amplitudes, convergence, and critical 
solutions. However the asymptotic dynamics of the bubble front in RMI appears quite 
different from that in RTI. In the RM family, the fastest solution corresponds to a bubble 
with a flattened shape: 


Ci = C a,rm, = 0 v = va,rm = 3/Akt (3.5) 

Higher order corrections for the solution (3.5) are reasonably small. The stability anal- 
ysis shows that solutions with a finite curvature C ~ 1/A are unstable, while flattened 
bubbles with ( = (a.rm are stable. Therefore, the physically significant solution in the 
RM family is the solution (3.5). 

The foregoing analysis can be applied for 3D flows with other symmetries. In either 
RTI or RMI, the nonlinear dynamics of 3D highly symmetric flows (Abarzhi 1998, 2002)) 
coincide except for the difference in the normalization factor k. A nearly isotropic shape 
of the bubble is the reason of this universality. The 3D and 2D results are similar quali- 
tatively. In main order, in RTI v a ,rt,3d/va,rt,2D ~ \/3 and Ca,rt,3D Ka,rt.2D ~ 3/4, 
while in RMI v a ,rm,3d/va,rm,2D ~ 2 and ( A ,rm,3d/Ca,rm,2D ~ 1, similarly to Abarzhi 
(1998) and Abarzhi (2002) for A = 1. For 3D low-symmetric flows, the asymptotic analy- 
sis shows a tendency of 3D bubbles to conserve isotropy in the plane, and a discontinuity 
of the 3D-2D dimensional crossover. 


4. Discussion 

Based on the foregoing results, we expect the following dynamics of the bubble front 
in the Rayleigh-Taylor and Richtmyer-Meshkov instabilities for fluids with finite density 
differences in the case of a small initial perturbation. In RTI, the bubble curvature Ci 
and velocity v grow as ~ exp (t/r R T) in the linear regime, t/r R T 1, and reach finite 
values Ci « C a,rt and v « v a ^ R t asymptotically for t/r R t ^ 1- The parameters (a,rt 
and v Ai rt depend strongly on the Atwood number. However, the value v\ RT /( A RT is a 
parameter universal for all A. This universality suggests that the bubble front evolution 
in RTI is quite complicated and cannot be approximated by the motion of a spherical 
bubble, in contrast to suggestions of the drag models of Oron et al. (2001) and Goncharov 
(2002). For A < 1 the bubbles are slower and less curved compared to the case of A = 1; 
for A « 0 the bubbles flatten and their velocity approaches zero. 

In RMI, the bubble curvature and velocity change as Ci ~ —kt/r R M and v(t) — Vq ~ 
—Vo t/r R M in the linear regime, t/r R M 1; for t ~ t R m the curvature Ci reaches an 
extreme value, dependent on the initial conditions; as t/r R M 1 the bubble flattens, 
Ci ~ C a,rm and decelerates v ~ v AR m- For A < 1 the bubbles move faster than those 
for A = 1, and for A < 1 the bubbles are flat. The flattening of the bubble front is a 
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Atwood number, A 



Figure 1. Dependence of the bubble curvature on the Atwood number A for 3D highly sym- 
metric flows in the Rayleigh-Taylor and Richtmyer-Meshkov instabilities; k is the wavevector, 
C,l = —Ak/8 is the Layzer-type solution, = —k / 8 is given by the drag models, (a,rt and 
C a,rm are given by the non-local multiple harmonic solutions in RTI and RMI respectively. 


distinct feature of RMI universal for all A. It follows from the fact that RM bubbles 
decelerate. The equation (3.5) suggests that for for fluids with very similar densities, 
A ~ 0, the bubble velocity has a much faster time-dependence, as t a with —1 < a < 0, 
in a qualitative agreement with experiments of Jacobs & Sheeley (1996). 

Figures 1, 2 and 3 compare our multiple harmonic non-local solutions with the models 
of Oron et al. (2001) and Goncharov (2002), and with the Layzer-type solution, which 
agrees with Sharp (1984). For 0 < A < 1 the asymptotic dynamics is different in RTI 
and RMI: the bubble velocity reaches a constant value in RTI and decreases with time 
in RMI. Therefore the distribution of pressure around the bubble is distinct in RTI and 
RMI. This should lead to a different shape for the bubble front. Our theory adequately 
describes the fact that in RTI the curvature C, a,rt has a strong dependence on the Atwood 
number, while in RMI Ca.rm = 0 for all A , Figure 1. The models of Sharp (1984), Oron 
et al. (2001), and Goncharov (2002) do not predict any difference between the shape of 
the RT and RM bubbles. In RTI, for A = 1 the values of Ca.rt = Cl = C d- However, 
for finite A the difference among the values Ca,rt , Cl and Cd is significant, Fig. 1; in 
the limit A — > 0, Ca,rt/Cd — > 0, while Ca,rt/Cl — > oo. We emphasize that for A = 0 
the RT instability does not develop and the bubble curvature should remain zero for all 
t, in agreement with (3.4). In RMI, the bubble in (3.5) flattens asymptotically for all A, 
while the models of Oron et al. (2001) and Goncharov (2002) suggest an A-independent 
finite value of the bubble curvature Cd = Cl\ a _ v Fig. 1. 

We conclude that in either RTI or RMI the bubble curvature is sensitive diagnostic 
parameter, which tracks the conservation of mass in the flow. The value A|£i| defines 
how flat or narrow the bubble is for a given lengths scale A. This dimensionless shape 
parameter is related to the bubble velocity and determines the flow drag. The value of 
the drag force is still a subject for controversy in the chaotic RTI and RMI (Dimonte 
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velocity, v/(g/k ) 1 12 



Figure 2. Dependence of bubble velocity v on the Atwood number A for 3D highly symmetric 
flows in the Rayleigh- Tay lor instability! h is the wavevector and g is gravity, vl^rt — \jAg/k 

is the Layzer-type solution, vd,rt = \f2Ag/(l + A)k is given by the drag models, and va,rt 
is given by the RT multiple harmonic non-local solution. 


2000). Our theory resolves this issue for highly nonlinear Rayleigh-Taylor and Richtmyer- 
Meshkov instabilities. The foregoing analysis suggests that the empiric models of Sharp 
(1984), Oron et al. (2001) and Goncharov (2002) may not be applicable in RTI and RMI. 

Our results are in good agreement with available data (Kucherenko et al. 2000; He 
et al. 1999; Gardner et al. 1988; Jacobs & Sheeley 1996; Holmes et al. 1995; Schneider 
et al. 1998; Holmes et al. 1999; Dimonte 2000; Volkov et al. 2001). However, in most of 
existing experiments, the measurement of the bubble curvature requires an improvement 
of diagnostics of the interface dynamics. The major issues to check by future observations 
are the following: Does the curvature of the RT bubble have a strong dependence on the 
Atwood number or reach an A-independent value? Does the curvature of the RM bubble 
vanish asymptotically or approach a finite value? 

For the bubble velocity, in the case of A ~ 1 the experiments and simulations have 
been in a reasonable agreement with the dependencies Vl,rt and vd,rt in RTI, and 
with the dependence Vd,rm in RMI. On the other hand, as evinced in the foregoing, 
the single-mode solutions Vl,rt and vd,rt as well as vl,rm and Vd,rm violate the 
conservation laws. This apparent paradox is easily explained. Figure 2 shows that in RTI 
for finite A and for £i ~ 1/A, the bubble velocity is relatively insensitive to details of the 
interface dynamics, and the quantitative distinction among the values of va.rt , vd.rt 
and VR'RT is 10 — 20%. In RMI, the velocities Vr^rm and Vr>,rm differ significantly from 
va,rm, Figure 3. The difference is however hard to distinguish in experiments. Since 
the velocity decays as 1/t, the bubble displacement, ~ Clog(t/rRM)/k, is comparable 
with an experimental error, and the coefficient C is in this dependence is impossible to 
evaluate. In contrast, the bubble curvature is a reliable diagnostic parameter as discussed 
in the foregoing. 

Our theory describes the principal influence of the density ratio on the large-scale 
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velocity, vAkt 



Figure 3. Dependence of the bubble velocity v on the Atwood number A for 3D highly 
symmetric flows in the Richtmyer-Meshkov instability; k is the wavevector, t is time, 
vl,rm = (2 — A 2 )/Akt is the Layzer-type solution, vd,rm = 2/(1 + A)kt is given by the 
drag models, va,rm = 3 /Akt is given by the RM multiple harmonic non-local solution. 


nonlinear dynamics of the RT and RM bubbles. The analysis is based on the assumptions 
that the flow dynamics is governed by a dominant mode, the transfers of energy to smaller 
or larger scales are not extensive, and the vorticity does not change the time-dependence 
of the large-scale coherent motion. If these conditions are broken (for example, for fluids 
with similar densities, A ss 0), the potential approximation may not give a correct time- 
dependence for the asymptotic bubble motion. We address these issues in the future. 


5. Conclusion 

We have found a multiple harmonic solution for a complete system of conservation laws 
describing the large-scale coherent dynamics in RTI/RMI in general three-dimensional 
case. The theory yields new dependencies of the bubble curvature and velocity on the 
density ratio, formulates the universal properties of the nonlinear dynamics in RTI and 
RMI, reveals an important difference between the dynamics of RT and RM bubbles, and 
shows the significance of mass conservation for the buoyancy-drag balance. 
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